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LATERAL STABILITY AND CONTROL DERIVATIVES OF A 
JET FIGHTER AIRPLANE EXTRACTED FROM FLIGHT. TEST DATA 
BY UTILIZING MAXIMUM LIKELIHOOD ESTIMATION 

By Russell V. Parrish and George G. Steinmetz 
Langley Research Center 

SUMMARY 

A method of parameter extraction for stability and control derivatives of aircraft 
from flight test data, implementing maximum likelihood estimation, has been developed, 
and successfully applied to actual lateral flight test data from a modern sophisticated jet 
fighter. This application demonstrates the important role played by the analyst in com- 
bining engineering judgment and estimator statistics to yield meaningful results. During 
the analysis, the problems of uniqueness of the extracted set of parameters and of longi- 
tudinal coupling effects were encountered and resolved. The results for all flight runs 
are presented in tabular form and as time history comparisons between the estimated 
states and the actual flight test data. 

INTRODUCTION 

A method of parameter extraction for stability and control derivatives of aircraft 
from flight test data has been developed at the Langley Research Center (ref. 1). This 
method, utilizing maximum likelihood estimation, has been applied to actual longitudinal 
flight test data from a modern sophisticated jet fighter airplane (ref. 2) to establish the 
merits of the estimation technique and its computer implementation. 

In the present study, the application of the method to actual lateral flight test data 
from the same airplane has also been used to establish the merits of the estimation tech- 
nique and its computer implementation by extracting, from the flight data, a set of stabil- 
ity and control derivatives that are well defined in terms of their standard deviations. 
During the analysis, the problems of uniqueness of the extracted set of parameters and of 
longitudinal coupling effects were encountered and resolved. The results presented 
demonstrate that the technique provides sufficient information to identify the uniqueness 
problem, if one exists, in terms of parameter correlations. The results also demon- 
strate that sufficient excitation of the aircraft will yield a unique set of derivatives and 
that incomplete modeling will be indicated by a poor fit of the data. 


The flight test runs utilized in this study are lateral responses generated by rud- 
der and/or aileron deflections in the neighborhood of ±10° and ±20°, respectively. The 
changes in angle of sideslip, roll angle, rolling velocity, yawing velocity, and lateral 
acceleration are typically ±10°, ±30°, and ±40° per second, ±15° per second, and ±0.2g, 
respectively. The parameters extracted were the standard linear body -axis lateral 
stability and control derivatives, with additional nonlinear derivatives dependent upon 
angle of attack. These nonlinear derivatives were found to be necessary due to strong 
longitudinal motion present in some of the flight test runs. 

SYMBOLS 


Measurements and calculations were made in the U.S. Customary Units. They are 
presented herein in the International System of Units (SI) with the equivalent values in 
the U.S. Customary Units given parenthetically. 


l Y,cg 


lateral acceleration at center of gravity, g units 


a. 


Y,I 


lateral acceleration at accelerometer location, g units 


wing span, meters (ft) 


mean aerodynamic chord, meters (ft) 


rolling-moment coefficient 


C *r 


damping-in-roll derivative, 


9C 


l 


i§) 


, per radian 


dC, 

C; = v per radian 

r #) 

ac 7 

C7„ effective -dihedral derivative, — per radian 

c /3 9/3 


9C, 

Ci. = — - per radian 
*6 a 96 a 


dC } 

Ci = — - per radian 
*6 r 96 r 
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rolling-moment coefficient at /3 = /3,p, 6 r = 6 r = T 

^T* r,T* a,T 


(Ci \ Cj at ii = oi— , per radian 

V L r) ant Lr T 


(%) a % at “ = <V per radian 

C n yawing -moment coefficient 

3C n 

Cn P = zm per radian 

\2V/ 

3C r 


C n „ = 


' n 


nr ~7^\ 

d \2v) 


per radian 


9C r 


: n p static directional -stability derivative, per radian 


, = ^n 

" n 6 a 96 a 

3C„ 


per radian 


% r = *r 

per 

(Ca K\T 

,6 a,T 

(S)„ T 

Cn P 

K)« T 

Cn ^a 


yawing -moment coefficient at /3 = f}^, 6 r = 5 r T’ 6 a = 5 a T 


C n at a = a—, per radian 


side -force coefficient 


8C, 


\2VJ 


dC-, 


per radian 


Cv = , , - per radian 


\2V> 


3 


C v = per radian 

Y 0 3/3 


9C 


Cvr = — - per radian 

96 


C Yk = 


9C-, 

96! 


per radian 


( C y)a « * side -force coefficient at /3 = 0 T , 6 r = 6 r t’ 6a = 6 a T 

' //3 T’ r,T’ a,T 

(%)« T % at ° - "T 

g acceleration due to gravity, meters/second^ (ft/sec^) 


ix 


aircraft moment of inertia about the body X-axis, kilogram -meters* 
(slug-ft^) 


L XZ 


lY 


iz 


K C, 


Kr 


Kr 


Kr 


Kr 


% 




' n 6 c 


L P 


product of inertia of aircraft referred to body X- and Z-axis, 
kilogram -meters^ (slug-ft^) 

p 

aircraft moment of inertia about the body Y-axis, kilogram -meters 6 
(slug-ft^) 

p 

aircraft moment of inertia about the body Z-axis, kilogram -meters 
(slug-ft^) 

p 

slope of linear variation of with a, per radian 

o 

slope of linear variation of with a, per radian 

slope of linear variation of C 0 p with a , per radian^ 

slope of linear variation of C n - with a, per radian^ 

°a 

p 

slope of linear variation of Cyp with a, per radian 6 


m 


mass of fueled airplane, kilograms (slugs) 
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p 

q 

r 

S 

u 

V 
v 
w 
Xy 

Y y 

z y 

a 

01 rp 

P 

/?T 

6 a 

V 


rolling angular velocity, radians/second 
pitching angular velocity, radians/second 
yawing angular velocity, radians/second 
wing area, meters^ (ft^) 

velocity along longitudinal body axis, meters/second (ft/sec) 

true airspeed, meters/second (ft/sec) 

velocity along lateral body axis, meters/second (ft/sec) 

velocity along vertical body axis, meters/second (ft/sec) 

accelerometer offset coordinate from center of gravity along longitudinal 
body axis, meters (ft) 

accelerometer offset coordinate from center of gravity along lateral 
body axis, meters (ft) 

accelerometer offset coordinate from center of gravity along 
vertical body axis, meters (ft) 

angle of attack, radians 

trim angle of attack, radians 

sideslip angle, radians 

trim sideslip angle, radians 

aileron deflection angle (positive when right aileron is deflected down), 
radians 

•aileron deflection angle at trim, radians 

rudder deflection angle (positive when trailing edge is deflected 
to the right), radians 
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6 r -p rudder deflection angle at trim, radians 

r fit error 

y,X arbitrary parameters 

6 pitch angle, radians 

c t > roll angle, radians 

p mass density of air, kilograms/meter^ (slugs/ft^) 

A dot over a variable indicates the time derivative of that variable. 

FLIGHT TESTS 

The flight test data were provided by the U.S. Naval Air Test Center at Patuxent 
River, Maryland. The flight tests were conducted by Navy test pilots as part of an inves- 
tigation with a McDonnell Douglas F-4 airplane. Five different lateral response runs 
were made: three during one flight test of the airplane and two during a second flight 
test. The first three runs were made at an altitude of approximately 6096 m (20 000 ft) 
at Mach numbers of about 0.6, 0.7, and 0.8, respectively. Control inputs for these runs 
were rudder only, rudder and aileron, and rudder only, respectively. The other two runs 
were made at an altitude of approximately 11 277.6 m (37 000 ft) at Mach numbers of 
about 0.9 and 0.8, respectively. Control inputs for these runs were rudder only and rud- 
der and aileron, respectively. The stability augmentation system (SAS) was deactivated 
in order to provide full response for all the test runs. 

For each of the test runs, the airplane was trimmed by the pilot at the desired alti- 
tude and Mach number and held for a short period. Then the control input or inputs were 
applied. No attempt was made to null any longitudinal motions. Roll and pitch angles as 
well as Mach number, pressure altitude, rudder deflection, aileron deflections, and cali- 
brated airspeed were recorded every tenth of a second. True airspeed was determined 
from figure 1 of reference 3 using Mach number, pressure altitude, and temperature 
from flight tests and resolved through angle -of -sideslip measurements to yield lateral 
velocity. 

Lateral displacement of the control stick in the F-4 airplane produces a combination 
of aileron and spoiler deflections. The aileron deflection is limited from 0° to 30° down- 
ward and from 0° to 1° upward. The spoiler being located on the upper surface of the 
wing has no downward deflection and is limited to upward deflections between 0° and 43°. 
In the flight records only the aileron deflections were recorded. Aileron-deflection data 
were used in the following manner to yield a single control input, which reflects a spoiler 


6 



effect. The assumption was made that a negative reading for either the right or left 
aileron was the indication of an aileron input. It was further assumed that the spoiler 
effect on the opposite side of the negative aileron deflection was equivalent to a positive 
aileron deflection of the same magnitude. Hence, by doubling the magnitude of the nega- 
tive aileron deflection and applying the sign convention of a right aileron to the magnitude, 
a single right aileron input, which is effectively the total aileron input, could be used in 
the equations. It should be noted that the aileron coefficients ^Cy g , C , and C nga J 

extracted by this program reflect the effect of both aileron and spoiler. Since these con- 
trol surfaces are physically linked, it is impossible to uniquely determine the coefficient • 
of each aileron and spoiler without additional information. 


Instrumentation consisted of rate gyros located slightly forward and at foot level of 
the pilot for measuring pitching, rolling, and yawing velocities; accelerometers located 
in the left wheel well for measuring lateral and normal accelerations; and vanes on a nose 
boom for measuring angle of attack and angle of sideslip. (See fig. 1.) No documentation 
was available from the Navy as to the accuracy of the instrumentation, although the method 
of parameter extraction (ref. 1) typically yielded the following signal-to-noise amplitude 
ratios (the noise amplitude was the 2 -sigma level): 


Lateral velocity 
Rolling velocity 
Yawing velocity 
Roll angle 

Lateral acceleration 


18 decibels 
24 decibels 
20 decibels 
22 decibels 
8 decibels 


AIRCRAFT MATHEMATICAL MODEL 


The equations of motion used by the computer program (ref. 1) were modified con- 
tinually during the analysis. However, three basic models evolved. The first model 
consisted of mainly lateral motion, the second model contained longitudinal coupling, and 
the third model contained longitudinal coupling and nonlinear lateral derivatives. The 
nonlinear derivatives Kq v , Kq^, Kq 7 , Kc n , and KQ nj _ permit variations with 


'l 


n 6c 


angle of attack in those particular derivatives that exhibit such dependence in wind-tunnel 
results (ref. 4). The three models can be obtained from the following equations: 


1 pV 


V = g cos 9 sin <p + pw - ru + — 

Cl 


( Cy p)a T +Kc Y p(“-“ T ) 



^T’ 6 r,T’ 6 a,T 


C Yp(P ~ ^t) + C Y 


rb 
r 2V 


' 6 *. - 6 , 


■,T) +C Y6 a ( 6 a 


6 ~ - 6 
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= p + tan 0(Q sin ^ + r C0S ^ 


4> * p 


(3 = sin_1 (?) 


V = \|u 2 + v2 + 


a 


= tan 




- i-(v + ru - P w " ® 

a Y,cg ' g 


cos e sin <P) 


a Y,I = ^» c e 


in uv - (p 2 + r 2 )Yy + fe* 

+ gUpq + ^ x y ' p 


- WZ^ 
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The first model, mainly lateral motion, can be obtained from the basic equations 
by requiring all longitudinal variables (u, w, 9, a, and q) to be constants and the non- 
linear derivatives Krv , K r , Kp, , Kp , and Kr- to be zero. The second 

Y P fcj3 W r ’ c n p ’ '-' n 6 a 

model, containing longitudinal coupling, can be obtained by using the longitudinal flight 
data as inputs to the equations in the same manner as rudder and aileron deflections are 
used. The third model, containing both coupling and nonlinear derivatives, is obtained 
from the second model by not restricting the nonlinear derivatives to zero. 

The nonlinear derivatives Kp , Kp, , Kp, , Kp„ , and Kp , as well as 

^Y p ’ 4 r ’ '-' n p c n 6a ’ 

the longitudinal coupling terms, were discovered to be necessary in order to fit the flight 
test data, as is demonstrated in the next section. Also demonstrated in the next section 
is the problem of uniqueness mentioned in references 1 and 2. 

RESULTS 

The conditions for the five flight test runs are listed in table I. The analysis of 
these runs involved two major problems: uniqueness and longitudinal effects. The 
results are presented in a manner to illustrate how each problem was encountered and 
then resolved. 


Uniqueness Problem 

Before presenting the uniqueness problem as encountered in this study, it would be 
well to describe the problem and the means of detecting its presence. The problem 
itself can be best described as follows: Given a set of parameters that minimize the fit 
error between measured and computed variables, does another set of parameters exist 
that will yield the same fit error? If the answer is yes, a uniqueness problem exists. 
Detection of the problem is facilitated by the use of the covariance matrix provided by 
the maximum likelihood estimation technique. Minor manipulation of this matrix, as 
described in reference 1, yields pairwise parameter correlation coefficients which esti- 
mate the degree of linear dependence between two parameters. Figure 2 illustrates the 
existence of a uniqueness problem due to linear correlation between arbitrary parameters 
X and y. Values of X and y that lie on the line of dependence yield the same fit 
error. However, it should be emphasized that two parameters may exhibit high correla- 
tion without indicating a uniqueness problem. Thus, it is necessary for the analyst to 
test any parameters with significant correlation coefficients to determine whether a 
uniqueness problem is present. The test is simply to determine whether the fit error 
changes as the parameters vary along the line of dependence. The procedure for carry- 
ing out the test is to assign to one of the correlated parameters several values in the 
range of interest and then extract the other parameter’s values; this determines the line 


9 



of dependence. In figure 2, the fit error r does not change as A and y vary along 
the line of dependence. Thus, a uniqueness problem is present. If the fit error did 
change, both parameters would be identified by the estimation technique at the point of 
minimum fit error and no uniqueness problem would exist, although the parameters 
would still be correlated. In this hypothetical illustration, the correlation between A 
and y is perfectly linear and will cause divergence of the estimation technique when an 
attempt is made to extract both parameters. However, in the use of real data, the pres- 
ence of noise usually prevents perfect linear correlation, and thus divergence. 

Figure 3 presents the model responses generated by the estimates of the stability 
derivatives of test run 1 and the respective flight test data, using the first model with all 
longitudinal variables fixed as constants (average values obtained from the flight data for 
each variable). (Note that symbols in figure 3 and subsequent machine plots presenting 
model responses and respective flight test data are not the standard symbols defined in 
the Symbols section.) Table II presents the estimates of the derivatives obtained, and 
table III presents a form of the covariance matrix for these estimates. Diagonal ele- 
ments of this matrix are the standard deviations of the estimates, and the off-diagonal 
terms are correlation coefficients. As denoted by the asterisks of table IE, Cy^, Cy p , 
and Cy r ; C^, C Zp , and C^; C nj3 , C np , and C nr ; and and C n ^ all have sig- 

nificant correlation. Investigation of these parameters revealed the existence of a unique- 
ness problem. 

A major cause of uniqueness problems is generally admitted to be insufficient 
excitation of the aircraft (for example, ref. 5). Test run 1 had rudder deflections only. 
Test run 2 contained both rudder and aileron deflections, and the model responses gener- 
ated by the derivative estimates for this test run and the respective flight test data are 
shown in figure 4. Again the longitudinal variables were fixed as constants during the 
extraction process. Table IV presents the estimates of the derivatives obtained, and 
table V presents the modified covariance matrix. As pointed out in section 7.8.3 of ref- 
erence 5, the likelihood of obtaining a unique set of derivatives is increased when both a 
rudder input and an aileron input are used to excite the airframe, as is evidenced by the 
lack of correlation exhibited in table V. 

Longitudinal Coupling Effects 

Examination of figure 4 (test run 2 responses) reveals poor fits for all the lateral 
variables; these poor fits indicate a possibly incomplete model. The longitudinal data 
for test run 2 are presented in figure 5 and indicate a substantial amount of longitudinal 
motion. Use of the longitudinal data as input, together with the modeling of angle of 
attack dependence of some of the derivatives, resulted in the extraction of a new set of 
derivatives for test run 2. Figure 6 presents the model responses generated by this set 
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of derivatives and table VI contains the derivatives and their standard deviations. No 
significant correlation was present and, thus, a unique set of derivatives has been 
extracted. It should be noted that a lack of confidence exists for all the Cy deriva- 
tives with the exception of Cy^, due to the large standard deviations of the estimates, 
as is the case with some of the nonlinear derivatives. 

Solution of the Uniqueness Problem 

The uniqueness problem of test run 1 was resolved by fixing the values of the non- 
linear derivatives and » C^, and Cn r at the values obtained in test run 2 

(the wind-tunnel results presented in ref. 4 show these derivatives to be fairly insensi- 
tive to Mach number variations in this flight regime) and extracting the remaining deriv- 
atives. This same procedure was used to solve the uniqueness problem of test run 3, 
which also had a rudder-only input. The model responses generated by the final esti- 
mates of the derivatives for test run 1 and test run 3 are shown with the respective flight 
data in figures 7 and 8, respectively. The values of the derivatives and their standard 
deviations are presented in table VII for test run 1 and table VIII for test run 3. 

Test run 4 had essentially a rudder-only input, whereas test run 5 had both rudder 
and aileron inputs. Again, the results of test run 5 were used to solve the uniqueness 
problem of test run 4. The model responses generated by the final derivative estimates 
of test run 4 and the respective flight data are shown in figure 9, and the estimates with 
the standard deviations are presented in table IX. The results of test run 5 are presented 
in figure 10 and table X. 

The total results of the analysis are summarized in figures 11 to 13, which illus- 
trate the variation of the extracted derivatives with Mach number, altitude, and angle of 
attack. The results shown in figure 13 are presented with the intercept values located 
at the trim angle of attack (symbol location) and the slope of the lines determined by the 
nonlinear derivatives. 


CONCLUDING REMARKS 

It is believed that the importance of the analyst, exercising engineering judgment 
tempered with estimator statistics, has been aptly demonstrated by the results of this 
study in recognizing and resolving the problems of uniqueness and longitudinal coupling 
effects. Thus, the extraction technique and its computer implementation have been shown 
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to provide the means for identifying both modeling and uniqueness problems and to yield 
a unique set of derivatives from actual lateral flight test data, provided the flight data 
contain sufficient information. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., August 4, 1972. 
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TABLE I.- FLIGHT TEST CONDITIONS 


Test 

run 

Altitude 

Mach 

number 

Center of 
gravity, 
% c 

Input 

m 

ft 

1 

6 096 

20 000 

0.6 

32.19 

Rudder 

2 

6 096 

20 000 

.7 

31.85 

Rudder and aileron 

3 

6 096 

20 000 

.8 

31.49 

Rudder 

4 

11 277.6 

37 000 

.9 

29.18 

Rudder 

5 

11 277.6 

37 000 

.8 

29.11 

Rudder and aileron 


TABLE II.- DERIVATIVE ESTIMATES OF TEST RUN 1 OBTAINED WITH MODEL 1 
Cy^ -0.392 


C Yp • • 
Cy r . . 

C Y 6r • 


C *r 


1.97 

3.75 

-0.0487 

-0.0938 

-0.355 

-0.230 


° n /3 

c n p 

Cn r 


0.00221 

0.120 

0.162 

-0.0664 

0.0462 
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TABLE IV.- DERIVATIVE ESTIMATES OF TEST RUN 2 OBTAINED WITH MODEL 1 


Cy^ -0.303 

Cv 0.858 

X P 

Cy r 2.18 

Cyg -0.00341 

Cvc - 0.0252 

Cl p -0.0452 

C 7 -0.126 

L P 

C lr 0.277 

C -0.0121 

C ; -0.118 

5 a 

% 0.105 

C n 0.0670 

n P 

C nr ■ -0.266 

C n6r 0.0519 
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TABLE VI.- FINAL DERIVATIVE ESTIMATES OF TEST RUN 2 


OBTAINED WITH MODEL 3 FOR a T = 0° 


Estimate Standard deviation 

C Y -0.341 0.0142 

r 

(C Yp ) a 0.086 0.271 

K Cyp 2.3- 5.60 

C Y r 1.27 0.301 

Cyr -0.035 0.0163 

°r 

C Y . 0.036 0.0180 

°a 

(Ci^ a -0.0382 0.000996 

K r , -0.73 0.0524 

*/3 

C, -0.363 0.00574 

-0.085 0.0175 

Kq 0.60 0.681 

Ci 5 -0.0163 0.000632 

C 7r -0.0350 0.000429 

C n/3 0.101 0.000791 

(Cn^ 0.0990 0.00801 

V P/a T 

K P -0.90 0.152 

u n p 

C nr -0.229 0.0164 

C nK 0.0527 0.000617 

“Oj. 

(C n6 \ -0.00074 0.000614 


Kp 0.04 0.0139 

C "«a 
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TABLE VII.- FINAL DERIVATIVE ESTIMATES OF TEST RUN 1 
OBTAINED WITH MODEL 3 FOR a T = 0° 

Estimate Standard deviation 


C Y0\ * • ' 

( Cy p)c t 

K Cy V 

Y p 

C Y r 

CY 6 r • 

(%)„ T 

K %' ' 

% 

( CZr )o! T ' 

K <V • * 

C ^r * * * 

% • • • • 

( Cn p)« T 

*C n 

n P 

Cn r 

Cn 6 r 

Fixed from test run 2. 


-0.295 

0.086 

2.3 

1.6 

-0.0240 

-0.0388 


-0.73 

-0.363 

-0.132 

0.60 

-0.0160 

0.0870 

0.0823 

-0.90 

-0.228 

0.0510 


0.0238 


* 

0.513 

0.0232 

0.000879 

* 

* 

0.0206 

* 

0.00132 

0.000393 

0.00291 


* 

0.000521 
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TABLE Vin.- FINAL DERIVATIVE ESTIMATES OF TEST RUN 3 
OBTAINED WITH MODEL 3 FOR a T = 0° 


C Y jS * 

(%)«. 

K % 

c Yr • 

Cy 6 r 

(%)«, 

K % 
% ' 

(%) ai 

C lS r ■ 

Cn (3 ' ' 

(SL 


C n 

u r 

C n 6r 

* Fixed from test run 2. 


Estimate Standard deviation 


-0.293 

0.0138 

0.086 

* 

2.3 

* 

1.20 

1.09 

-0.0513 

0.0164 

-0.0486 

0.000643 

-0.73 

* 

-0.363 

* 

-0.0318 

0.0143 

0.60 

* 

-0.0164 

0.000648 

0.0987 

0.000491 

0.0974 

0.00423 

-0.90 

* 

-0.228 

* 

0.0483 

0.000618 
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TABLE IX.- FINAL DERIVATIVE ESTIMATES OF TEST RUN 4 
OBTAINED WITH MODEL 3 FOR a T = 0° 



Estimate 

Standard deviation 

% 

-0.269 

0.0164 

( Cy p)c t 

0.090 

* 

K % 

2.3 

* 

C Y r 

1.03 

0.578 

C Y 6r * 

-0.0767 

0.0155 

CY «a 

0.0602 

* 

( C ^ T 

-0.0549 

0.000979 

% 

-0.367 

* 

C; 

-0.282 

* 

( C : *r)a T 

0.170 

0.0244 

K 'V r 

0.80 

* 

Ci 6 r 

-0.00363 

0.000978 

C ^ a 

-0.0244 

* 

Cn i3 

0.107 

0.000533 

(° n p)o T 

0.125 

0.00369 

K s 

-0.90 

* 

Cn r 


* 

Cn 6 r 

0.0542 

0.000506 

( Cn 6 a ) 

' a 'Qf«p 

-0.00529 

* 

Kr_ 


* 


n 6 a 


* Fixed from test run 5. 
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TABLE X.- FINAL DERIVATIVE ESTIMATES OF TEST RUN 5 


OBTAINED WITH MODEL 3 FOR a T = 0° 


Estimate 

-0.383 

Standard devi; 
0.0434 

0.090 

0.181 

2.3 

* 

1.01 

1.21 

-0.0473 

0.0342 

0.0602 

0.00931 

-0.0577 

0.00752 

-0.367 

0.0875 

-0.282 

0.00531 

0.0362 

0.0410 

0.80 

0.824 

-0.00360 

0.00145 

-0.0244 

0.000682 

0.0965 

0.00107 

0.0987 

0.00526 

-0.90 

0.493 

-0.158 

0.0351 

0.0525 

0.000907 

-0.00529 

0.000342 

0.04 

* 


Fixed from test run 2. 
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Figure 2.- Hypothetical illustration of correlation of two parameters. 
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Figure 3.- Model responses generated by derivative estimates of test run 1 
with longitudinal variables constant. 
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Figure 8.- Model responses generated by final derivative estimates of test run 3 

with longitudinal data as input. 
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Figure 8.- Concluded. 
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Figure 10.- Continued. 
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Figure 12.- Variation of control derivatives with Mach number. 
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